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ABSTRACT 

We report the measurement of the three-point correlation function (3PCF) 
of galaxies for the Las Campanas Redshift Survey (LCRS). We have not only 
measured the 3PCF in redshift space but also developed a method to measure 
the projected 3PCF which has simple relations to the real space 3PCF. Both 
quantities have been measured as a function of triangle size and shape with 
only a fractional uncertainty in each individual bin. Various tests derived from 
mock catalogs have been carried out to assure that the measurement is stable 
and that the errors are estimated reliably. Our results indicate that the 3PCFs 
both in redshift space and in real space have small but significant deviations 
from the well-known hierarchical form. The 3PCF in redshift space can be 
fitted by Qred{s,u,v) = 0.5 ■ io[o-2+o-i{7tt)']^' for 0.8 < Si2 < 8h-^Mpc and 
S31 < 16h~^Mpc, and the projected 3PCF by Qproj{rp,u,v) = 0.7r~^-^ for 
0.2 < rpi2 < 3 /i~^Mpc and rp3i < 6 /i~^Mpc (s and are in unit of /i~^Mpc), 
though a systematic weak increase of Qproj{rp,u,v) with v at Vp > 1 h~^Mpc is 
noted. The real-space Q{r,u,v) for 0.2 < ri2 < 3 /i~^Mpc and r^i < Gh^^Mpc 
can be well described by half the mean 3PCF predicted by a CDM model with 
Qoh = 0.2. 

The general dependence of the 3PCF on triangle shape and size is in 
qualitative agreement with the CDM cosmogonic models. Quantitatively the 
3PCF of the models may depend on the biasing parameter and the shape of 
the power spectrum, in addition to other model parameters. Taking our result 
together with the constraints imposed by the two-point correlation function and 
the pairwise velocity dispersion of galaxies also obtained from the LCRS, we find 
that we have difficulties to produce a simple model that meets all constraints 
perfectly. Among the CDM models considered, a flat model with Q = 0.2 meets 
the 2PCF and PVD constraints, but gives higher values for the 3PCF than 
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observed. This may indicate that more sophisticated bias models or a more 
sophisticated combination of model parameters must be considered. 



Subject headings: galaxies: clustering - galaxies: distances and redshifts - 
large-scale structure of Universe - cosmology: theory - dark matter 



Introduction 



Correlation functions are very powerful statistics to describe the large scale structures 
in the Universe ( [Peebles 1980, hereafter P80). The lowest order, the two-point correlation 
function (2PCF) ^(r), has been widely used to measure the clustering strength of galaxies 
and to confront models of cosmic structure formation. Quite a number of large galaxy 
catalogs, both angular and redshift, have been used to determine the 2PCF, and this 
statistic has now been established quite well ( |Jing, Mo, fc Borner| 1998, hereafter JMB98; 
Lin et aL| 1998 in preparation; [Lin et aE] 1997; [Tucker et al.[ 1997; [Ratclitfe et al.| 1997; 
Baugh[ 1996; [Hermit et al.[ 1996; [Strauss fc Willie^ 1995 for a review before 1995). This 
statistic has produced several constraints on theoretical models already despite the fact 
that there are many ingredients to a specific model which can be optimally adapted to 
the properties of the galaxy sample. The cosmological parameters, like the initial power 
spectrum of the DM component and the bias, i.e. the difference in the clustering of galaxies 
and DM particles, can all be adjusted to some extent. 

The three-point correlation function (3PCF) C('^i25 ^"235 ''^31) is a further statistic useful 
in characterizing the clustering of galaxies (P80). Its measurement can give additional 
constraints for cosmogonic models. The determination of the 3PCF was pioneered by 
Peebles and his coworkers in the seventies. Based on their careful analysis of the Lick and 
Zwicky angular catalogs of galaxies they propose a so-called "hierarchical" form 



C(n2,'^23,'^3i) = Q i{ri2)i{r2z) + ^(^^23)^(^31) +^(^^31)^(^^12) 



(1) 



with the constant Q ~ 1.29 ± 0.2. This form is valid for scales r < 3/i~^Mpc (P80). The 
analysis of the ESO-Uppsala catalog of galaxies ( [Lauberts| 1982) by [Jing, Mo, fc Borner 
(1991) supports this result. The 3PCF was also examined for the CfA, AAT and KOSS 
redshift samples of galaxies ( [Peebles[ 1981; Pean et aT] 1983; [Efstathiou fc Jedrzejewsk 



1984; Pale-Sutton et ah] 1989). Because all these redshift samples are small (with < 2000 
galaxies), these authors were not able to examine the validity of the hierarchical form in 
redshift space. Instead they just forced a fit of the hierarchical form and obtained the value 
of Q. The Q value of redshift samples obtained in this way is around 0.6 ( [Ffstathiou fc 
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Jedrzejewski 1984), much smaller than the value advocated by Peebles and his coworkers. 
The difference may partially be attributed to the redshift distortion effect which reduces 
the Q value ( |lVlatsubara| 1994). The skewness analysis of the 1.2 Jy IRAS survey of galaxies 



has given a similar Q value ( |Bouchet et al.| 1993). 

The hierarchical form (Eq. |l]) is purely empirical. There is no solid theoretical argument 
supporting this form. In contrast, the second-order perturbation theory predicts that Q 
depends on the shape of the triangle and on the slope of the linear power spectrum 



1984) in the linear regime. For a CDM-like power spectrum with a slope which changes 
with the scale, Q then varies with the size and shape of a triangle ( |Jing fc Borner| 1997). 
Even in the strongly non-linear regime where the hierarchical form was established, the 
CDM models do not seem to obey this form as demonstrated by [Matsubara fc Sut"o| (1994) 
based on N-body simulations. Recently, |Yano fc Gouda| (1997) have re-examined, based on 
the BBGKY equation, the stable clustering (strongly non-linear) problem and pointed out 
that the hierarchical form holds if the clustering is stable. The question is if the condition 
of stable clustering can be achieved in the real Universe (|Jain| 1997). 

The 3PCF of galaxies carries much useful information which is important for 
cosmogonic models. The theories based on CDM models predict that the 3PCF of galaxies 
depends on the shape of the linear power spectrum (|Fry| 1984; |Jing fc Borner| 1997) and 



the galaxy biasing relative to the underlying mass ( Pavis et al| 1985; [Gaztanaga fc Freiman 



1994; |Mo, Jing, fc White] 1997; [Matarrese, Verde, fc Heavens| 1997; [Catelan et aT| 1998). It 



might also be sensitive to a possible non-Gaussianity of the initial density fluctuation ( Pry 



Scherrei] 1994). Furthermore the 3PCF must be determined accurately if one wants to 



use the cosmic virial theorem to obtain the mean density of the Universe. 

We have carried out a detailed analysis of the Las Campanas Redshift Survey 
( [Shectman et al.| 1996), and in this paper we report the measurement of the 3PCF of 
galaxies in this survey. We have not only measured the 3PCF in redshift space but also 
developed a method to measure the projected 3PCF which has simple relations to the real 
space 3PCF (§3.1). Our methods are checked very carefully with the help of mock catalogs 
generated from N-body simulations, and the physical meaning of these two quantities is 
also investigated (§3.2). Our statistical results are compared with previous work, with 
the emphasis on a critical examination of the hypothesis of the hierarchical form (§3.4). 
As we will see, the hierarchical form does not seem to be a good prescription even in the 
strong clustering regime. We will present a new fitting formula in §3.3. Implications for 
cosmogonic models are discussed in §4. 
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2. Observational sample and mock catalogs 

The sample used for our analysis is the Las Campanas Redshift Survey ( |Shectman| 
et al.| 1996; hereafter LCRS). This is the largest redshift survey, which is now publicly 



available. Our sample consists of all galaxies with recession velocities between 10,000 and 
45,000 kms~^ and with absolute magnitudes (in the LCRS hybrid R band) between —18.0 
and —23.0. There are 19558 galaxies in this sample, of which 9480 are in the three north 
slices and the rest in the three south slices. The survey is a well-calibrated sample of 
galaxies, ideally suited for statistical studies of large-scale structure. All known systematic 
effects in the survey are well quantified and documented (Shectman et al. 1996; Lin et al. 
1996), and so most can be corrected easily in statistical analyses. The only exception is 
the 'fiber collision' limitation which prevents two galaxies in one ~ 1.5 x 1.5 deg^ field from 
being observed when they are closer than 55" on the sky, because it is impossible to put 
fibers on both objects simultaneously. Here we will use extensively mock catalogs generated 
from N-body simulations to quantify this effect. 

The real-space 2PCF and the Pairwise Velocity Dispersion (PVD) have been determined 
for the LCRS by JMB98. The redshift-space 2PCF and power spectrum for this sample 
were presented by [lucker et al.| (1997) and [Lin et al] (1997) respectively. All these studies 



have shown that the LCRS is large enough to accurately measure these low order statistical 
quantities. In particular, JMB98 have carried out a detailed comparison between the 
observed 2PCF and PVD and the predictions of currently favoured CDM cosmogonies. 
They have used a large set of mock samples to adequately compare models and observations. 
The construction of mock catalogues from the simulations, i.e. photometric catalogues 
subject to the same selection effects as the real observations are a very important aspect 
of their analysis, because only in this way could the statistical significance of the results 
be asserted. Three spatially fiat models have been considered in JMB98, with [Qq, Aq, 
r,a8) = (0.2,0.8,0.2,l.), (0.3,0.7,0.2,1.), and (1.0, 0.0, 0.5, 0.62), where Qq is the density 
parameter, Aq is the cosmological constant, F = Q^h and ag are the shape parameter and 
normalization of the CDM power spectrum (Bardeen et al. 1986). All of the models give 



a steeper 2PCF, and a higher PVD on small scales than the data. Thus unless galaxies 
are biased with respect to the mass with a scale-dependent bias, all these models can be 
ruled out. Unfortunately physical models for a density or a (not so wanted, but perhaps 
unavoidable) velocity bias are not on firm grounds. Therefore in JMB98 a simple, but 
plausible phenomenological model for the bias has been suggested. To suppress the number 
of pairs in the DM distribution at small separations, it is assumed that the number of 
galaxies per unit dark matter mass N/M is smaller in massive halos than in less massive 
ones. If a behaviour such as N/M oc with a = —0.08 is used for clusters of mass M^, 
the predictions of some CDM models are consistent with the observational results. The 
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best agreement was achieved for the flat fig = 0.2 modeL 

We will use 10 mock catalogs of this model to test our statistical methods and 
quantify the 'fiber collision' effect. Since this model has reproduced the LCRS 2PCF 
and PVD, we believe these mock catalogs are very suitable for this purpose. We shall 
also use these mock samples for model testing, as an example to illustrate the power of 
the three-point correlation function in discriminating between models which have similar 
two-point correlations. Since the model is a typical CDM model, we will generalize the 
discussion to other CDM models. 



The three-point correlation function (3PCF) C('"i2, ''"23, '"si) is defined, through the joint 
probability dPi23 of finding one object simultaneously in each of the three volume elements 
dri, dv2 and dr^ at positions ri, r2 and respectively, as follows (P80): 



dPi2z = n(ri)n(r2)n(r3) [1 + ^{ru) + ^(r23) + ^(rsi) + ((^12, ^23, r3i)]dri(ir2(ir3 (2) 



where rij = |rj — rj|, and n{ri) is the mean density of galaxies at r^. This definition can be 
applied straightforwardly to redshift surveys of galaxies to measure the 3PCF C(si2, S23, S31) 
of galaxies in redshift space (at this point we neglect the anisotropy induced by the redshift 
distortion which will be considered later). Here and below we use r to denote the real space 
and s the redshift space. 

The 3PCF of galaxies can be measured from the counts of different triplets (P80). 
For this purpose, a sample of randomly distributed points, which has exactly the same 
boundaries and the same observational selection effects as the real survey, is generated. 
Four types of distinct triplets with triangles in the range (si2 ± l/2Asi2, S23 ± I/2AS23, 
and S31 ± I/2AS31) are counted: the count DDD{si2, S23, S31) of triplets formed by three 
galaxies; the count DDR{si2, S23, S31) of triplets formed by two galaxies and one random 
point; the count DRR{si2, S23, S31) of triplets formed by one galaxy and two random points; 
the count RRR^Su, S23, S31) of triplets formed by three random points. Following the 
definition [eq(H)], we shall use the following estimator 



3. The three point correlation function 



3.1. 



Definitions and statistical methods 



C(si2,S23,S3l) 



27RRR\si2,S23,S3l) X DDD{si2,S23,S3l) 
DRR^{si2, S23,S3l) 

9RRR{si2, S23, S31) X DDR{si2, S23, S31) 



+ 2 



(3) 



DRR'^{Si2,S23,S3i) 
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to measure the 3PCF of the galaxies in redshift space. The above formula is slightly 
different from the estimator used by Groth & Peebles (1977). Here we have extended the 
argument of piamilton| (1993) for the 2PCF to the case of the 3PCF. The coefficients 27 and 



9 are due to the fact that only distinct triplets are counted in this paper. Since the early 
work of Peebles and coworkers (P80) indicates that the 3PCF of galaxies is approximately 
hierarchical, it is convenient to express the 3PCF in a normalized form Qredisu, S23, S31): 

^ / X _ C{Si2, S23, S31) 

^red[Sl2, S23, S31 — — r—- ^ , 777 TTTT \C? ^ " ^'^ 

USi2)US23) + US23)US3i) + ^(s3l)USi2) 

It is also convenient to use the variables introduced by Peebles (P80) to describe the 
shape of the triangles formed by the galaxy triplets. For a triangle with the three sides 
■S12 < S23 < S31, s, u, and v are defined as: 

s = S12, u = — , v = . (5) 

S12 S12 

Clearly, u and v characterize the shape and s the size of a triangle. We take equal 
logarithmic bins for s and u with the bin intervals Algs = Algw = 0.2, and equal linear 
bins for v with Av = 0.2. For our analysis, we take the following ranges for s, u and v: 
0.63 < s < 10/i^^Mpc (6 bins); 1 < m < 4 (3 bins); and < < 1 (5 bins). 

A sample of 25,000 random points is first generated. The counts RRR are less than 
~ 5 for small triangles (s < 1 /I'^Mpc). In order to suppress the fiuctuation induced by the 
random samples, we have recalculated the counts RRR for S31 < 4 /I'^Mpc by generating 
a random sample 10 times larger, which ensure that the counts RRR are at least ~ 300 
of the interested triangle configurations. We scaled these counts to 25,000 random points 
and also use these counts to get DRR on the small scales since RRR/ DRR is constant. 
However, it is not trivial to search triplets for so many points. We have generalized the 
ordinary linked-list technique of P^M simulations (Pockney fc Eastwood] 1980) to spherical 
coordinates to count the triplets. The linked-list cells are specified by the spherical 
coordinates, i.e. the right accession a, the declination 5 and the distance s. With this 
short-range searching technique, we can avoid the triplets out of the range specified thus 
making counting triplets very efficient. 

The 3PCF in redshift space Qred{s,u,v) depends both on the real space distribution of 
galaxies and on their peculiar motions. Although this information contained in Qred{s, u, v) 
is also useful for the study of the large scale structures (see §4), it is apparent that 
Qred{s,u,v) IS different from Q{r,u,v) in real space. In analogy with the analysis for the 
two-point correlation function, we have determined the projected three-point correlation 
function H(rpi2, rp23, Tp3i). We define the redshift space three-point correlation function 
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Cz{rpi2, rp23, rp3i, 7ri2, Trig) through: 

dPl23 ^ ^(Sl)n(s2)n(s3)[l +^^(rpi2,7ri2) + 6('^p23, 7r23) +^z{rp3UT^3l) 

+Cz{rpi2, rp23, r-psi, nu, Tii3)]dsids2ds3 (6) 

where rfPf23 the joint probabihty of finding one object simultaneously in each of 
the three volume elements dsi, ds2 and dss at positions Si, S2 and S3; is the 

redshift space two-point correlation function; Vpij and TTjj are the separations of objects 
i and j perpendicular to and along the line-of-sight respectively. The projected 3PCF 
n(rpi2,rp23,Tp3i) is then defined as: 

n(rpi2, rp23, rp3i) = j (2(^2, ^^p23, rp3i, ttu, 'Ki3)dni2d'K23 (7) 

Because the total amount of triplets along the line-of-sight is not distorted by the peculiar 
motions, the projected 3PCF n(rpi2, rp23, ^^^31) is related to the 3PCF in real space 
((^^12,^23,^31) : 

n(rpi2, rp23, rp3i) = j ((^'^^12 + ^12, \/^p23 + y23, + {yi2 + y23Y)dyi2dy23 ■ (8) 

Similarly as for C(si2, 523, •S31), we measure G(V2, '^p23, '^p3i, 7ri2, 7713) by counting 
the numbers of triplets -D-D-D(rpi2, rp23, ^^^31, 7ri2, 7ri3), DRR{rpi2,rp23,rp3i, '^12,^^13), 
RDD{rpi2,rp23,rp3i,ni2,'Ki3) and RRR{rpi2,rp23,rp3i, 7:12, 7:13) formed by galaxies and/or 
random points with the projected separations rpu, rp23, and rp3i and radial separations ir^ 
and 7123. We will use rp, u and v: 

r p23 1^p31 — fp23 /^x 
Tp = rpi2, U = , V = . (9) 

^pl2 '"pl2 

to quantify a triangle with rpi2 < rp23 < rp3i on the projected plane. Equal logarithmic bins 
of intervals Algr^ = Algii = 0.2 are taken for rp and u, and equal hnear bins of Av — 0.2 
for V. The same ranges of u and v are used as for ^(s,u,v), but is from 0.128 /i~^Mpc 
to Ah^^Mpc (7 bins). The radial separations 7ri2 and 1123 are from — 25 /i~^Mpc to 
25 h-^Mpc with a bin size of 1 h ^Mpc. The projected 3PCF is estimated by summing up 
(^(rp, 7r^2, """23) at different radial bins (vr|2,7r23): 

U(rp, u,v)^Yl Cz{rp, u, v, nl^, 71^3) A7ri2A7r^3 (10) 

and normalized as 

P' ' w{rpi2)w{rp23)+w{rp23)w{rp3x)+w{rp3i)w{rpi2)' 
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where w{rp) is the projected two-point correlation function ( Pavis fc Peebles 1983; JMB98) 

w{rj,) = Y,Urp,7^')A7r' (12) 

i 

An interesting property of the projected 3PCF is that if the three-point correlation function 
is of the hierarchical form, the normalized function Qproj{fp,u,v) is not only a constant 
but also equal to Q. Therefore the measurement of Qproj{rp, u,v) can be used to test the 
hierarchical form which was proposed mainly based on the analysis of angular catalogs. 



3.2. N-body tests of the statistical methods 

To test the reliability of our statistical analysis and to demonstrate the effects of 
the redshift distortion, the projection, and the fiber collisions, we make use of the full 
simulation and the mock catalogs. Because the mock samples are cluster-(under)weighted, 
we have applied the same weighting to the full simulation to achieve a proper comparison. 
To calculate the quantities in redshift space for the full simulation, we assume that the 
third axis is along the line-of-sight. 

In Fig. (H) we compare Qred{s,u,v) estimated with our statistical method from the 
mock samples with the true value. The latter is determined from the full simulation using 
the method of [Jing fc Bornei) (1997). On the scales from 1 h~^Mpc to 10/i~^Mpc, the two 



estimated quantities agree fairly well, indicating that the LCRS can yield an unbiased 
estimate within the error bars of the 3PCF in redshift space. The test is important, 
considering the fact that the LCRS is essentially two-dimensional with one dimension in 
the direction of the line-of-sight. 

Fig. (H) shows the projected Qproj{rp,u,v) estimated from the mock samples with the 
method described above. The true Qproji^p, u, v) can be calculated from the real space 3PCF 
({r,u,v) through Eq. (|]). We determine ({r,u,v) for the full simulation using the method 
of |Jing fc Borner| (1997) and calculate the integral of Eq.(^ by linearly interpolating the 
estimated ({r,u,v). The two estimated quantities agree very well within the error bars, 
indicating that our method can give a correct estimate of the projected 3PCF. The real 
space 3PCF Q{r,u,v) is also shown in the figure by the thick lines. It decreases with the 
scale r, as noted previously ([Matsubara fc Suto| 1994; |Jing fc Borner| 1997). The reason for 
the decrease is due to the fact that the slope of the power spectrum is more negative on 
smaller scales (see Jing, in preparation, for a detailed discussion). The consequence is that 



r„, M, f IS 



ojy pi 



due to the averaging of Q{r,u,v) on scales r > Vp [Eq.(|D], the projected Qpr, 
also a decreasing function of but smaller than Q{r,u,v) for r = Vp (compare thick and 
thin lines in the figure). Another interesting point is that Q{r,u,v) is much higher than its 
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counterpart in redshift space Qred{s,u,v) on scales < 10/i~^Mpc. This result is well-known 
since the redshift distortion smears out the dense clusters and thus reduces the 3PCF on 
small scales (e.g. [Matsubara fc Su"to| 1994; |Matsubara| 1994). 

We have also tested for the fiber collision effect of the LCRS sample as we did for the 
2PCF and PVD in JMB98. While both the 2PCF [^(s) and w{rp)] and the 3PCF [({s,u,v) 
and n(rp, u, v)] show some small dependence on this effect, it cancels out completely, when 
we divide these two quantities by one another to form the normalized 3PCF Qred{s,u,v) 
and Qproj 

{rp,u,v). Therefore, the fiber collisions of the LCRS have little effect on the 
normalized functions Qred{s,u,v) and Qprojijp^u^v). 

In summary, these tests have convinced us that our method is suitable for giving a 
stable measurement of the 3PCF from the LCRS with reliable error estimates. 



3.3. The statistical results of the LCRS 

We present our results of the 3PCF in redshift space Qred{s,u,v) and of the projected 
3PCF Qproj{rp,u,v) in Figures (|^) and @ respectively for the Las Campanas Redshift 
Survey. The errors of the Q-values are the bootstrap errors which are estimated with the 
approximate formula of |Mo, Jing, fc Bornei| (1992). As we can see from Fig. (|^), the 3PCF 
obtained from redshift space is not changing very much with s or u; it increases somewhat 
with V. For small v, Qred is approximately constant with a value of ~ 0.5, but it increases 
up to ~ 1 when v ^ 1. Compared with the 3PCF in redshift space, the projected one 
Qproj{rp,u,v) [Figure (^] shows quite similar dependences on triangle shape (i.e. u and v), 
but quite a different dependence on the triangle size (i.e. s or Vp), though the errors of the 
projected 3PCF are larger. Its value decreases with Vp from about 1.2 at Vp = 0.2 h^^Mpc 
to 0.5 at Tp about 2/i"^Mpc. Both this decrease with growing rp and the weak increase 
with V are in contrast to the hierarchical assumption. If the three-point correlation function 
in real space were hierarchical, equation (P) shows that the projected one would also be 
hierarchical and equal to Q. This behavior, however, is qualitatively in agreement with the 
CDM model predictions, and we will discuss this point in §4. 

The hierarchical form (eq.l) does not seem to provide an adequate description of our 
results of the 3PCF in redshift space Qred{.s,u,v) or of the projected 3PCF Qproj{i^p,u,v). 
We have looked for fitting formulae for both quantities. Our results can be fitted quite well 
by Qred{s,u,v) = 0.5 ■ lo'°'^"^°'^^^^ and Qproj{rp,u,v) = 0.7r~°-^ (s and rp are in unit 
of h~^Mpc), which are the smooth lines on Figure (§) and the thick lines on Figure (|^). 
We have neglected the weak systematic dependence on v of Qproji^p, u, v) at rp > 1 h^^Mpc 
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in fitting this quantity since the dependence is not statistically significant on the scales we 
probed. We do not intend to give an error estimate for the coefficients in the formulae, 
since the errors of the Q-values in individual bins are likely non-Gaussian distributed and 
correlated among different bins. These formulae are intended to give a simple, but for most 
purposes accurate, approximation to our statistical results. 

Since the real space 3PCF ({r,u,v) possibly depends on r, u and w in a complicated 
way, the inversion of eq.(|]) to get ({r,u,v) from the projected Il{rp,u,v) is certainly 
unstable. We have noted however that the projected Qproj{rp,u,v) can be modeled very 
well by half the value of the mock projected Qpf^j{rp,u,v) [the thin lines in Fig. (|)], which 
means that Q{r,u,v) = 0.5Q'^'^"^{r,u,v), where Q'^'^"^{r,u,v) is the real-space 3PCF of the 
CDM model, is a good approximation to the real-space 3PCF of LCRS galaxies on scales 
< 3 /i^^Mpc. These Q{r, u, v) can be easily read out from the (5'^'^™(r, u, v) in Fig. (|^). 



3.4. Discussion 

The high quality of the LCRS survey, by its CCD photometry, its complete redshift 
information and its large size, has enabled us to give a reliable determination of the 3PCF. 
It is the first time that the three-point correlation functions both in redshift space and in 
real (projected) space can be measured as a function of the triangle size and shape, with 
only a fractional error in each individual bin, for a wide range of triangle configurations. 
Although our statistical results in real space on scales < 1 /i~^Mpc are not far from the 
hierarchical prediction (1) with Q = 1.3 ±0.2 ( proth fc Peebles] 1977; GP77), the systematic 
changes of Qproj{rp,u,v) with the triangle configurations clearly point to a more elaborate 
model for Q{r,u,v). One such model was already proposed in §3.3. Although our result 
of Qred{s,u,v) generally agrees with the previous studies based on much smaller redshift 
samples (see §1), our measurement has much better accuracy. 

Recently there are concerns about the reliability of the high-order correlation functions 
derived from photographic-plates based galaxy catalogs. The skewness of the galaxies 
derived from the Automatic Plate Machine (APM) angular catalog ( |Maddox et al.| 1990) and 
from the Edinburgh/Durham Southern Galaxy Catalogue (EDSGC, peydon-Dumbleton,| 
[Collins, fc IVlacGillivrayl 1989) is significantly different even on small scales < 1 h^^Mpc 
( |Gaztahaga| 1994; [Szapudi, IVleiksin, fc INicEo]] 1996), despite the fact that both catalogs 



are constructed from the same UK Illa-J Schmidt photographic plates and the latter is 
just a sub-sample (~ 1/4) of the former. The difference seems to arise from their different 
methods to digitalize the plates ( Szapudi fc Gaztanaga 1998), and at this point it is rather 



difficult to judge which method is superior. It is also worth noting that the 3PCF results 
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of GP77 are somewhat sensitive to the corrections they have apphed to the Lick catalog. 
Fortunately the LCRS survey does not suffer from these uncertainties since its photometric 
catalog was constructed from the CCD drift scans. Furthermore in our analysis we did not 
have to apply any additional corrections except for those well-documented by the survey 
team. 

It is interesting to compare our results with the skewness Ss{R) determined from the 
large angular catalogs, in particular since the two studies based on the APM and EDSCG 
catalogs have yielded rather discrepant results. The skewness is related to the 3PCF 
through an integral as follows: 

^'<^) - m 

C{R) = tt:t c?rirfr2(ir3C(ri2,r23,r3i) , 

V'^ Jspherofi 
1 



an) = TT^ / dr.dr^aru), (13) 

V J sphere R 

where V = ^R^. It needs full information of Q{r,u,v) to calculate the skewness [eq. ([T3|)], 
and we use the model proposed in §3.3. The skewness 5*3 (i?) for the LCRS survey is then 
about 4.5 a.t R = 0.2 h~^Mpc and about 3.5 at i? = 1 h^^Mpc, which seem in agreement 
with the results of [Szapudi, IVleiksin, fc lNichol| (1996) based on the EDSGC catalog but 



significantly higher than the APM results of |Gaztahaga| (1994) on the scales < 1 h ^Mpc. 



The cosmic virial theorem (CVT) has been widely used to measure the mean density 
of the universe. If the 3PCF is hierarchical, the CVT can be expressed in its simplified 
form relating the density parameter Qq, ^(r), Q, and the PVD cri2(r). Since our results 
show that the 3PCF is not hierarchical, this relation becomes much more complicated 
and it is necessary to work out the integration over C('"i2 5 ''"23, ''"31) which might depend on 
C(^i2, ''^23, ''^31) on very small scales r ~ 0. Since C('^i2? ''"23? ''^31) is a decreasing function of the 
triangle size, previous studies which usually used the Q value at ~ 1 h~^Mpc might have 
overestimated the mean density. The size of galaxies, which are usually treated as point 
sources in the CVT application, may also be important in the estimate of the mean density, 
especially near r ~ ( Peebles] 1976; ^uto fc Jing| 1997). 



We would like to remark here that the significant difference between Qred and Qproj 
at small scales comes from peculiar motions of galaxies. This gives another possibility to 
estimate the velocity dispersion of galaxies ( |lVlatsubara| 1994). 



4. A case for model testing 



In this section we compare the 3PCF of the LCRS with model predictions. |Jing fc 
Borner| (1997) and Jing (in preparation) have recently studied the 3PCF Q{r,u,v) for a set 
of CDM models based on second-order perturbation theory and N-boby simulations. We 
found that for fixed u and v, Q{r, u, v) is a decreasing function of the size r. In the strongly 
non-linear regime {C,{r) ^ 1, Q{r,u,v) shows a very weak dependence on u and v. In the 
weakly non-linear and linear regimes, Q{r,u,v) increases significantly with v for fixed r 
and u. All these features are found in the projected Qproj{rp,u,v) of the LCRS galaxies. 
Therefore the statistical results found from the LCRS survey are all qualitatively consistent 
with the mass 3PCF based on N-body simulations of cosmological models. 

As an example to quantitatively test models with the 3PCF, we compare the 3PCFs 
of the LCRS galaxies with the results of the mock samples in Figs. (|^) and (^. Only in 
this way could the redshift distortion and the projection effects be accounted for properly. 
From the figures we find that the qualitative features, i.e. the dependence on v for fixed 
s or Tp, and u, the decrease of Q with increasing values of s or Vp are reproduced quite 
well in the mock samples. The values of the data set, however, are consistently lower than 
the mean model predictions by a factor ~ 2. Since the Q- values in each bin of the ten 
mock samples are not Gaussian distributed (skewed to high values), it is not meaningful to 
use the standard deviation to quantify the statistical significance. Instead we pick up the 
lowest of the ten mock Qred or Qproj values in each bin and compare it with the observed 
results. The thick lines in Figs. (|^) and (P) correspond to these lowest values for u = 1.29, 
which should be compared with the open triangles. These lines are still higher than (in 
most bins) or at least as high as (in a few bins) the observational values, which indicates 
that the observational values are lower than this model's predictions at a confidence level 
> 90%. The underlying model for the mock sample is a CDM universe with Qq = 0.2 
and Ao = 0.8 and with clusters underweighted (see §2). We have computed the 3PCF 
for this universe without cluster weighting, and found the cluster weighting, like the fiber 
collision effect, does not change much the value of Q{r,u,v). Thus, even if this model fits 
the 2PCF and the PVD of the LCRS galaxies quite well, it seems not really adequate to 
describe the clustering of galaxies when the 3PCF is considered, unless the observed 3PCF 
is biased low by the cosmic variance [at < 10% probability] . This might indicate that the 
gravitational interaction alone is not sufficient to describe the clustering of galaxies, and 
physical processes of gas and radiation hydrodynamics connected with galaxy formation 
must be taken into account. A positive bias, i.e. a biasing parameter 6 > 1, can reduce the 
3PCF. But, perhaps this conclusion goes too far - in fact, a slightly higher shape parameter 
r will give a better fit, because Q becomes smaller if F is increased. It appears that a model 
with a new set of parameters is to be sought and the 3PCF determined here should provide 
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to such new models a test in addition to the 2PCF and the PVD. 



5. Conclusion 

The result is clear, and the conclusions are straightforward: We have succeeded in 
measuring the 3PCF from the LCRS. This is the first time that the three-point correlation 
function of galaxies has been measured accurately from a redshift survey. Both the 3PCF 
in redshift space and the projected 3PCF have been measured as a function of the triangle 
size and shape with only a fractional uncertainty in each individual bin. Various tests have 
been carried out to assure that the measurement is stable and that the errors are estimated 
reliably. Our results indicate that the SPCFs both in redshift space and in real space have 
small but significant deviations from the well-known hierarchical form. The 3PCF in the 
redshift space can be fitted by Qred{s,u,v) = 0.5 ■ lof°-^+°-^(^)']^' for 0.8 < Si2 < 8 /i^^Mpc 
and S31 < 16/i~^Mpc, and the projected 3PCF by Qproj{^piU,v) — 0.7r~°-^ for 
0.2 < rpi2 < 3/i~^Mpc and rp3i < 6/i~^Mpc (s and Vp are in unit of h'^Mpc). Although it 
might be not unique to get 3-D real space Q{r, u, v) from the measured projected function 
Qproj{fp-i u, v), we found that a half of the predicted Q{r, u, v) of the CDM model considered 
in this paper provides a good description of the LCRS data. 

The three-point correlation function gives an additional statistical tool to constrain 
cosmogonic models. The general dependence of the 3PCF on triangle shape and size is 
in qualitative agreement with the CDM cosmogonic models. Quantitatively the 3PCF of 
the models may depend on the biasing parameter and the shape of the power spectrum, 
in addition to other model parameters. Taking our result together with the constraints 
imposed by the two-point correlation function and the pairwise velocity dispersion of 
galaxies also obtained from the LCRS, we find that we have difficulties to produce a simple 
model that meets all constraints perfectly. Among the CDM models considered, the flat 
model with Q — 0.2 meets the 2PCF and PVD constraints, but gives higher values for the 
3PCF than observed. This may indicate that more sophisticated bias models or a more 
sophisticated combination of model parameters must be considered. 

We are grateful to Yasushi Suto for helpful discussions, and for the hospitality extended 
to us at the physics department of Tokyo university. G. B. thanks the Yamada foundation 
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Comparison of mock to full simulation 
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Fig. 1. — The normalized 3PCF in redshift space Qred{s, u, v) of the mock samples (symbols) 
and of the full simulation (lines). The error bars are the la standard deviation of the 
measurement for the 10 mock samples. For clarity, the error bars are plotted lor u — 2 only 
but those for the other two values of u are very similar. 
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Comparison of mock to full simulation 
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Fig. 2. — The normalized projected Qproj{rp,u,v) of the mock samples (symbols) and of 
the full simulation (thin lines). The latter is computed through eq. (^ from the 3-D 3PCF 
Q{r,u,v) of the full simulation which are also plotted with the thick lines. A comparison of 
the symbols and the thin lines tests the statistical methods, and a comparison of the thin 
lines and thick lines shows the projection effect. As in Fig. (1), the error bars are the la 
standard deviation of the measurement for the 10 mock samples and are plotted for m = 2 
only. 
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Redshift-space 3PCF of the LCRS Survey 
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Fig. 3. — The normalized 3PCF in redshift space Qred{s, u, v) of the LCRS survey (symbols). 
The errors are estimated by the bootstrap resampling method. For clarity, the error bars are 
plotted for u = 2 only but those for the other two values of u are very similar. The results 
are weU fit by Qred{s, u, v) — 0.5 • iQ^^-^'^'^-^^ifi^ 1*' (s is unit of /i~^Mpc) which are shown by 
the solid lines. 
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Projected 3PCF of the LCRS Survey 
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Fig. 4. — The normalized projected 3PCF Qproj{i^p,u,v) of the LCRS survey (symbols). 
The errors are estimated by the bootstrap resampling method. For clarity, the error bars 
are plotted for u = 2 only but those for the other two values of u are very similar. The 
results are well fit by Qproj{i^p, u, v) = 0.7(rp)^°'^ (r^ is unit of /i^^Mpc) which are shown by 
the thick sohd hues. The thin lines are a half ol the mean projected 3PCF of the the mock 
samples, which seem to fit the observational data very well. 
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Comparison of the LCRS to the Model 
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Fig. 5. — The normalized 3PCF in redshift space Qred{Sj u, v) of the LCRS survey (open 
symbols) compared with the mean value of the ten mock samples (solid symbols). The thick 
solid lines represent the lowest of the ten mock Qred{s,u,v) values in each u — 1.29 bin, 
which can be regarded as the lower model hmits at ~ 90% significance level. 
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Comparison of the LCRS to the Model 
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Fig. 6. — The normalized projected 3PCF of the LCRS survey (open symbols) compared 
with the mean value of the ten mock samples (solid symbols). The thick solid lines represent 
the lowest of the ten mock Qred{s, u, v) values in each u — 1.29 bin, which can be regarded 
as the lower model hmits at ~ 90% significance level. 



